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Vacuum or matter free space is the most widely 
studied electromagnetic medium, not only because 
it underlies the development of continuum electro- 
magnetic properties from microscopic principles [1, 
2], but also because of the significance of electro- 
magnetic communication devices in modern society 
[3]. The electromagnetic constitutive equations of 
vacuum are commonly stated in textbooks as 

D(r,t) = eoE(r,i), (1) 
B(r,t) =/xoH(r,t), (2) 

where eo = 8.854 x 10"-^^ F m~^ and fiQ = Air x 
10~^ H m~^ in SI units, whereas r and t indicate 
position and time. 

These equations presuppose either the absence 
of a gravitational field or that the observer is local. 
When a gravitational field is present, spacetime ap- 
pears curved — which is well-known [4]. One can 
still use the textbook versions of the Maxwell pos- 
tulates for gravitationally affected vacuum, but the 
constitutive relations are now [5, 6] 



T>{r,t)=eo2ir,t) •E(r,t) 

-Coir(r,t) xH(r,t), 



(3) 



B{r,t) = ^^o2{r,t) '-tliv,t) 

+ Coir(r,i) xE(r,i), (4) 

in lieu of (1) and (2). Here 7(r, t) is a real symmet- 
ric dyadic and r(r, t) is a vector with real- valued 
components, both related to the metric of space- 
time; whereas co = 

Just as isotropic dielectric-magnetic mediums 
provide material counterparts of (1) and (2) in the 
frequency domain [1, 7], the vast variety of com- 
plex materials [8, 9] — natural as well as artifi- 
cial — suggests that it is quite possible that (3) and 
(4) also have material counterparts. This thought 
inspired the present communication, wherein we 
present the derivation of the dyadic Green func- 
tion for frequency-domain electromagnetic fields in 
a homogeneous medium inspired by (3) and (4). 

With the assumption that all fields have an 
exp(— iwf) time dependence, with oj as the angular 
frequency, the constitutive relations of the chosen 
medium are 

D(r, (j) = eo 7(0;) • E(r, oj) 

-Co'T{cu)xn(r,Lo), (5) 
B(r,w) = /xo7(w) • H(r,w) 

+ Co^T{Lo)xF,{r,u;). (6) 



The coordinate system has becm chosen such that 
7(0;) is diagonal, and from now onwards the depen- 
dence on w is implicit. 

The frequency-domain Maxwell curl postulates 
in the chosen medium may be set down as 



V X E(r) =iu! Ho 2' H(r) + Cq^ T x E(r) 



(7) 



V X H(r) 



eo7-E(r)-c-irxH(r) 



(8) 



where J(r) is the source electric current density. 
Our objective is to find the dyadic Green functions 
G (r, s) and G (r, s) such that 

E(r) = iw/xo y" y" j a^{r,s) '3{s)d^s,{9) 

H(r) = J J J ajr,s)-3{s)dh, (10) 

with the integrations being carried out over the re- 
gion where the source electric current density is 
nonzero. 

To begin with, the substitution of (9) in (7) and 
comparison of the resulting expression with (10) 
yields 

ajr, s) = 2"' • (V X I - tko r X i) • £^(r, s) , 

(11) 

where fco = uj^/eojlo and / is the identity dyadic. 
Thus, an expression for only G (r, s) has to be 
found. 

For that purpose, following Lakhtakia & Welgl- 
hofer [10], we start by defining new fields and source 
current density as 

e(r) = E(r) exp(-ifco T • r) , (12) 
h(r) = H(r) exp(-ifco Ft), (13) 
j(r) = J(r) exp(-iA;or -r). (14) 

Hence, (7) and (8) respectively transform to 

V X e(r) = za;/io7 • h(r) , (15) 

V X h(r) = -ia;eo7 • e(r) -l-j(r). (16) 

Next, we make use of an affine transformation 
associated with the scaling of space as per [11] 



f = 7V2 . 



(17) 
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where ^^^"^ • ^^^"^ = 7 and we recall that 7 is a 
dyadic with real- valued elements. Let us define an- 
other set of fields and source current density as 

e(r) = 7I/2 . 6(71/2 . r) , (ig) 
h(r) = 7I/2 . h(7i/2 . r) , (19) 



j(r) = adj7 



,,1/2 



■ j(7 



1/2 . 



(20) 



where 'adj' stands for the adjoint. Then, (15) and 
(16) transform to 

V X e(r) = iujuo g h(r) , (21) 

V X h(r) = -iujeo g e(r) + j(r) , (22) 

where 

9 = ^\, (23) 
and I7I denotes the determinant of 7. 

Prom the foregoing equations, we obtain 
[(V X I) • (V X I) - klg^^ • V X e(r) 

= iu) Hog j{r) . (24) 
The solution of (24) is well-known as [12] 

e{r) ^ iuJiiog j j j g{r,s) • i{s) d^s , 

where 

VV \ exp(-<fco5|r - s|) 



(25) 



3(r,s)= / + 



4n\r — s| 



(26) 



In order to go back from (25) to (9), we have to 
invert the two transformations in reverse sequence: 
Substitution of (18) and (20) in (25) yields 



e(r) = iw/io (adj V ^ 

g{r,s) 'j{s)d^s 



(27) 



wherein 



£(r,s) = (L + 
exp 



-7 . VV 



ikogh~^^^ ' (r - s) 



47r|7-i/2 . (r - s) I 



(28) 



By substituting for e and j in (27) using (12) and 
(14), respectively, we find 

E(r) = iw/io exp(iA;or • r) ^adj7^/^ j 
£(r,s) • J(s) 

exp(-iA;or • s)d^s). (29) 



Therefore, the dyadic Green function G^(r, s) emerges 
from (29) as 

g^(r,s) = exp[ifcor • (r- s)] (adj^'^^^ 



>^og — 



7 • VV 



exp 



*fco.9|7-'/' • (r-s)l 



47^1^-1/2 . (r - s) 



• (30) 



Equation (30) is the desired result. If 7 = / and 
r = 0, this expression reduces to the usual dyadic 
Green function for gravitationally unaffected vac- 
uum [12]. 
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